# plot histogram of interview density
hist10 <- hist(ESS10$date, breaks = 500,
               main="",
               xlab="Days from Taliban takeover",
               cex.lab = 1.5,
               cex.axis = 1.5,
               cex.main = 3)


# colored plot 
Treatment <- as.factor(ESS10$afg)

hist <- ggplot(ESS10, aes(x = date, fill = Treatment)) +
  geom_histogram(binwidth = 1) +
  theme_light(base_size = 25) +
  xlab("Days from Taliban takeover") +
  ylab("Interview Frequency")
hist
